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Abstract 

Peculiarities of ferromagnetic resonance response of conducting magnetic bi-layer films of nano- 
metric thicknesses excited by microstrip microwave transducers have been studied theoretically. 
Strong asymmetry of the response has been found. Depending on the order of layers with respect 
to the transducer either the first higher-order standing spin wave mode, or the fundamental mode 
shows the largest response. Film conductivity and lowered symmetry of microwave fields of such 
transducers are responsible for this behavior. Amplitude of which mode is larger also depends on 
the driving frequency. This effect is explained as shielding of the asymmetric transducer field by 
eddy currents in the films. This shielding remains very efficient for films with thicknesses well below 
the microwave skin depth. This effect may be useful for studying buried magnetic interfaces and 
should be accounted for in future development of broadband inductive ferromagnetic resonance 
methods. 
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I. INTRODUCTION 



Broadband inductive microwave setups are emerging as a common tool with which to 
measure dynamical properties of magnetic thin films and nano-structures jl,y, ^, Q]. Their 
main part is a section of a coplanar or a microstrip microwave transmission line ('"microwave 
transducer") on which top the sample under investigation sits. Absorption of microwave 
power by the transmission line loaded by the sample is measured by a network analyzer ^ 
or by a combination of a microwave generator and microwave passive components |^]. These 
setups are useful for accurately measuring spin wave excitations which are used as a probe 
for characterizing magnetic surfaces and buried interfaces. However, consideration needs to 
be made of the experimental setup when analyzing results . In particular, 

Schneider et al. [5| mention that measured response is dependent on the distance of the 
film from the coplanar surface due to a change in strip line impedance. In another work 
they report on using a floating ground plane above the film in order to enhance the response 
amplitude [l^. The broadband inductive technique was also shown to excite not only the 
fundamental ferromagnetic resonance mode, but also higher order spin wave resonances ^. 

In this paper we present theoretical results showing asymmetric broadband-FMR response 
of multilayered metallic films with total thicknesses much smaller than the microwave skin 
depth. This effect is one more feature which distinguishes this technique from the cavity 
ferromagnetic resonance (FMR). 

On the other hand, recently an interest has emerged to excitation of propagating spin 
waves in metallic magnetic films by narrow (several microns in width) microstrip microwave 



transducers 
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14j . This interest is largely due to the possibility of using spin waves 



to perform logic operations 
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3,Q, 



20|. If necessary, the expressions we obtain 



can be easily extended to describe excitation of propagating spin waves following suggestions 



in Ref . 



2l| and Ref. 



II. DYNAMIC EQUATIONS 



To construct our model we use the system of Maxwell equations with no electric-bias 



term as suggested previously 



V X h = ^ e 



(1) 
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V ■ h = -V • Anm 

V X e = -i^^(h + 47rm) 

Here m is the dynamic magnetization, e is an electric field, and uj is the eigenfrequency 
of spin waves. To obtain the last of these equations we assumed that all dynamic fields and 
the dynamic magnetization have the time dependence in the form exp{iujt). 

In Ref . jol starting with these equations we constructed a quasi-analytical theory allowing 
for inhomogeneous dynamic magnetization in the film plane in the form of spin waves with 
a set of in-plane wave numbers kx. This was important in that work, since due to the type 



of symmetry in t 
precession kx=Q 



le coplanar transducers, their microwave fields do not excite homogeneous 
24| . On the contrary, microstrip transducers show maximum of microwave 



absorption at kx=^ (see e.g. [25| and references therein) due to symmetry of the in-plane 
component of the dynamic magnetic field of the transducer hx- For microstrip transducers 
with a microstrip width much larger than the free spin-wave propagation path (about 100 
micron) in Permalloy (which is the metallic magnetic material with the lowest spin wave 
damping) one can consider the dynamic magnetization and all microwave fields being ho- 
mogeneous in the film plane (the (x,z) plane). Assuming that the microwave current of 
the microstrip and the static magnetic field H (which will later enter the magnetic torque 
equation) are along the axis these equations can be reduced to a system as follows: 

dhx/dy = -—^(^z (2) 

dejdy = i—{hx + 47rm^) 
c 

The last of these equations shows that the out-of-plane component of the dynamic mag- 
netic field is purely of magnetostatic nature and represents the dynamic demagnetizing field 
of niy. It remains present when the sample conductivity is set to zero. On the contrary, 
from the first two equations one sees that the in-plane component of the magnetic field is 
of curling nature. Indeed, one can reduce Eqs.(2a) and (2b) to a second-order differential 
equation: 

- i6^d'^hx/dy'^ + hx = -47rm^. (3) 



where S = \J c/ {Aiiauj) is the microwave skin depth. One sees that in the limit ma;=0, when 
the sources of the magnetostatic (potential) field vanish this equation remains valid. On 



the other hand, from Eqs.(l) one can find that hx vanishes for o" = 0. Thus this component 
is not present in the insulating samples, where all fields are known to be magnetostatic, it 
is entirely of curling nature, and is due to the eddy current jz induced in the conducting 
ferromagnetic films by the precessing magnetization and the external microwave magnetic 
field. 

In Ref. ^1 it has been shown how to obtain a particular solution of such inhomogeneous 
equation in terms of Green's functions. A general solution of the homogeneous equation has 
to be added to this solution to obtain the complete solution. This is not difficult to do, 



if magnetic layers composing the multi- layered film have the same conductivity [26|]. But 
if interfaces of magnetic metal layers with different values of conductivities are present the 
final analytic expressions become too cumbersome to be used in computations. 

Therefore, instead of solving Eq.(3) analytically we derive boundary conditions for this 
equation and solve it numerically together with Landau-Lifshitz-Guilbert equation. From 
the condition of continuity of and at the interface of two magnetic layers y = Li we 
find 

a^dh^^^/dy = a,dh^^/dy (4) 

where the indices 1 and 2 denote the layer numbers. The boundary conditions at the 
outer surfaces of the bilayer film are obtained by first allowing the outer space having finite 
conductivity (Tout- Then the magnetic field in the outer space is described by Eq.(3) with the 
zero right-hand part. Its solution is obtained requiring vanishing of the microwave magnetic 
field on both infinities. At one side from the film y < we have: = Aexp{^/i6outy), 
and at the other side from the film y > Li + L2 we have = Bexp{—^/i6outy)- Obviously 
the same boundary conditions (4) are valid at the boundary of a conducting magnetic film 
with a conducting outer space. Substituting the second of these expressions into Eqs.(4) one 
obtains a formula which relates the field inside the ferromagnetic layer to its derivative. 

1 CTout dhx ^ ^ ,_^ 

K = TT- at y = Li + L2, (5) 

ViSout ^2 oy 

where (T2 is the conductivity of the second layer of the bi-layer film. To derive the formula 
which relates to dh^/dy at the other film boundary we need to include the magnetic field 
induced by the microstrip. We model the microstrip line as a surface current density jz at 



the film surface y = 0. Then following the same procedure one obtains: 

1 out 9hx , V 

hx = ^ 3z at ?/ = 0. (6) 

In the limit aout « 0"i,cr2 Eqs.(5) and (6) reduce to: 

hx = -jo at y = 0, (7) 
K = at y = Li + L2. 

(In these expressions we use SI units.) The second of these formulas shows that the 
dynamic magnetic field at the film surface not facing the microstrip line is zero. This suggests 
that total back-reflection of the microwave magnetic field from the boundary between two 
media with a large difference in electric conductivity takes place. From the second of Eqs.(4) 
one now finds that the dynamic magnetic field outside the film y > Li + L2 vanishes. Thus 
a conducting magnetic layer with a thickness much smaller than 5 efficiently shields the 
microwave magnetic field. This is in a striking contract with insulating films: the microwave 
magnetic field of a transducer easily penetrates through magnetic insulators. It equals —jz/'^ 
everywhere for y > and jz/'i for y < 0. (see e.g. Eq.(32) in 251]). (To obtain this result 
from our expressions one has to set ai = (J2 = (Tout and then take the limit of vanishing 
conductivities.) The latter result allows one to consider the total field as a sum of the 
external field —jzf^ and the field of eddy currents. Then one finds that the field of the eddy 
currents should grow from —jz/'i at y = to jz/2 at y = Li + L2 ensuring no total dynamic 
magnetic field outside the film y < and y > Li + L2. 

The obtained boundary conditions (7) are quite asymmetric. This is a big contrast 
compared with the cavity resonance where a homogeneous external microwave magnetic 
field penetrates the sample through both surfaces and the shielding effect is not seen in the 
boundary conditions (Eq.(4.1) in 3]). (We note, that the conditions (7) are in agreement 
with the theory [22] and for the cavity resonances can be transformed into Eqns. (3.3) or 
(4.1) in the cited work.) 

In this work we numerically solve the linearized Landau-Lifshitz-Gilbert magnetic torque 
equation [27 1 

ium = (8) 
-7[(i^ma, + iyrriy + Mi^) 



where M is saturation magnetization of a layer, A is the exchange constant for the layer, 
He = H + iauj/'-f, a is the Gilbert magnetic damping constant, and i are the unit vectors 
along the coordinate axes. (The static magnetic field H and the equillibrium magnetization 
are along the z-axis.) The torque equation is solved together with Eqs.(2c),(3), the electro- 
magnetic boundary conditions Eqs.(7), and the exchange boundary conditions at the layer 
interface. The latter are cast in the form suggested in Provided there is no pinning 
of magnetization in both layers at the interface these exchange boundary conditions at the 
interface y = Li oi two layers 1 and 2 read 

9mW/5, + ^m«-^^m(^) = (9) 

where Au is the inter- layer exchange constant. For the out-of-plane dynamic magnetization 
components the inter-layer boundary conditions are the same. The exchange boundary 



conditions at the outer surfaces of the film are, as follows 



dnix/dy ± drrix = 0, drUy/dy = 0, (10) 

where d is the pinning parameter at the film surface. The positive sign is for the boundary 
y = and the negative one is for y = Li + L2. To obtain the numerical solution the 
system of differential equations is transformed into a system of finite-difference equations. 
The latter form a matrix-vector equation with a band matrix of coefficients. This linear 
algebraic system is solved using numerical methods of linear algebra. The way we incorporate 
boundary conditions between layers and at the outer surfaces of the film into the equations 
is shown in the Appendix. 

III. TRANSDUCER RESPONSE 

There are two ways to treat the transducer response. One of them is that of the effective 



microwave susceptibility 



]. It assumes that the total microwave energy absorbed 
by the material is lost due to magnetic losses by driven precession of magnetization in the 
area of localization of the microwave magnetic field of the transducer. Another method 
was established several decades ago when magnetostatic spin waves in monocrystalline fer- 
rimagnetic films of yttrium iron garnet (YIG) were found being promising for processing 



microwave signals [30|. It was suggested that linear impedance could be used to quanti- 
tatively characterize efficiency of excitation of propagating spin waves by microstrip (see 



e.g 



251, 131|, M, 1331, Cl and references therein), and coplanar [24| transducers. The linear 
impedance is calculated as a ratio of the Poynting vector of the flux of energy of microwave 
field through the transducer surface to the microwave current in the transducer 25(]. This 
approach remains valid also in the case when the transducer width is much larger than the 
spin wave free propagation path and the most of the irradiated energy is not carried away by 
spin waves but is lost due to relaxation processes within the reach of the microwave Oersted 
field of the transducer. In this case it naturally incorporates possibility of energy losses due 
to irradiation of propagating spin waves which leave the area of the transducer's microwave 
Oersted field and carry energy away. This additional loss mechanism for metallic rnagnetic 
films was recently discussed in Ref.j2| and experimentally studied in detail in 11, 12, 13|. 



Both the effective microwave susceptibility and the linear impedance can be related to the 
scattering coefficient 5*2 1 measured in the network- analyzer based broadband FMR j^. 

In this work we choose the complex linear impedance Zr of a microstrip loaded by the 
ferromagnetic film as a quantity describing the efficiency of microwave absorption. This is 
because in this way, if necessary, our theory can be easily extended to include the finite 
width of the transducer to describe effects of irradiation of propagating spin waves with 
nonvanishing wave numbers by narrow transducers which were used in recent experiments 
[ll , Q, [3] ■ Ref. 21 1 shows the way the irradiation of propagating waves can be treated. 

Our calculations are carried out following the suggestion in 25|| originally made for mag- 
netostatic spin waves in YIG films. We adopt this approach to metallic ferromagnetic films. 
Note, that the latter couple to the transducers much more weakly and take much less power 
from the transducer. For the in-plane homogeneous surface current and the in-plane homo- 
geneous microwave electric field the expression for Eq.(27) from 25| reduces to: 

e^iy = 0) 



(11) 



where w is the transducer width in the direction x. 

The obtained values of Zr are then transformed into the value of the scattering coefficient 
5*21. We start with the formula for the input impedance Zf of a section of a microstrip line 
loaded by a magnetic film. The film sample has a length / along z and sits on top of the 
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transducer. Following Eq.(25) in 25(] we have 



2^2^ zocosh{-ffl) + Z^sinh{{-ifl)) 
^ ZQsinhij fl) + ZcCosh{{jfl)) 

In our case of the broadband FMR zq is the characteristic impedance of the sections of the 
microstrip line not covered by the sample ("unloaded microstrip") and equals 50 Ohms. 
is the characteristic impedance of the section of the microstrip line loaded by the sample 
("loaded microstrip") 



Z, = ^{Z, + Z,)/Y^ (13) 

with Zq and Yq being the complex series resistance and the complex parallel conductance of 
the unloaded microstrip, and 7/ is the complex propagation constant of the loaded microstrip 



7; = ^(Zo + Z,)Y,. (14) 
The transmission matrix T of the loaded microstrip is obtained following Ref. 35l |: 



T= [T« ■T(') -T^^^)], (15) 

where T^^) and T^^) are the transmission matrices of junctions of the loaded microstrip with 
the unloaded microstrip. The former is for the junction at the front edge of the sample 
and the latter is for the rear edge. These matrices are defined via the complex reflection 
coefficient: 

r = ±|i-f . (16) 

where the positive sign is for the front edge and the negative sign is for the rear edge. The 
elements of these matrices are: 

TSl'^)=T?/) = (l-r)-\ (17) 

tS'/) = tS'^^ = r(i-r)-^ 

The transmission matrix for the loaded microstip between these two edges T^^-* has only 
diagonal elements: 

tSi^ = 1/t5? = e^/'. (18) 

The scattering parameter 5*21 of the whole loaded microstrip is I/T22 (Eq.(15)). By multi- 
plying the matrices in Eq.(15) one obtains: 

- 1 

521 = — — J. (19) 



The maErnitude of the hnear radiation impedance of magnetostatic spin waves in YIG is 

n 

of order of zq [2J] and one has to use Eqs. (12-14) as they are. But for nanometric metaUic 
films and wide microstrip transducers Zr << 50Q which is clearly seen in experiment as 



r |« 1 



36|. This allows considerable simplification of these formulas. From the 



condition << Zq to the first order in Z^j Zq from Eqs. (12-20) one obtains 

S21/521o = exp[-|^/], (20) 
where 5*210 is the scattering parameter of the transducer with no sample on its top. 



IV. DISCUSSION 

Results of numerical calculation using this formalism for the frequencies ^7/(277) 4, 7.5, and 
18 GHz are shown in Fig. 1. We consider a Cobalt (Co)-Permalloy (Py) bi-layer film with 
the parameters shown in the figure caption. The bilayer has a thick Py layer and a thin Co 



layer which role is to introduce dynamic magnetization pinning [37| of magnetization in the 
Py layer at the layers' interface while avoiding formation of additional resonances localized 
in the Co-layer. 

From the figure one sees that the absorption when the Co-layer faces the microstrip is a 
few times smaller than when the Permalloy layer faces it. One also sees that the amplitude 
of the first standing spin wave (the second peak from the right) is visible only for Co facing 
the transducer. It grows with frequency and becomes larger then the fundamental mode 
(the most right-hand peak) at higher frequencies. This result agrees well with experiment 



which will be published elsewhere 



36|. 



The next figure (Fig. 2) explains this strong asymmetry of the film response. Its upper 
panel (Fig. 2(a)) shows distribution of dynamic magnetization mx{y) across the bi-layer 
thickness when the surface current is applied from the side of the thin Cobalt layer. One 
observes partial dynamic pinning of magnetization at the Py-interface with Co which man- 
ifests itself as considerably inhomogeneous m2;(y)-distribution through the Py- layer for the 
fundamental mode with the minimum at the layer interface. The first higher-order mode 
of the stack represents a combination of the fundamental mode of the Co-layer and the 1st 
SSW of the Py-layer. For the surface current applied at y = Li + L2 one obtains the mirror 
image of this panel. 
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Figure 2(b) shows the /ia;(y)-dependence for both cases of film orientation with respect to 
the microstrip. First one sees that the distribution satisfies the boundary conditions Eq.(7), 
and in order to satisfy them the magnetic field needs to have a large negative gradient 
through the film. This field may be thought as a combination of the Oersted field of the 
microstrip and of the shielding field of the eddy current je in the film. The former is constant 
through the film and the space y > Li + L2 and equals —jo/2. The latter equals —jo/ 2 
at y < 0, grows through the film to reach +jo/2 at y — Li + L2, and remains equal to 
this value for y > Li + L2 to ensure no dynamic magnetic field behind the film. From this 
consideration on can infer that the direction of the eddy current is opposite to jg. If one 
looks at the plot of the electric field (Fig. 2(c)) one finds that the eddy current jg = ^^&z 
has a phase of tt and thus is indeed anti-aligncd to Jq. 

From Fig. 2(b) one also sees that hx{y) is practically a linear function. Indeed, the 
deviation from the linearity for the long solid line is less than 10 percent at the maximum 
of deviation. If one looks at Eq.(3) one finds that without the right-hand part this equation 
has the solution h% satisfying the boundary conditions (7) as follows: 

K^A sinh[Vi(5(y - Li - L2)] (21) 

(here we neglect the inter-layer boundary) . In the case case of a thin monolayer conducting 
film Li « 5, L2 — Eq.(22) reduces to a hnear dependence — A\/i5{y — Li), where 
A — 2jo/[exp{2\^SLi) — 1] jo/{ViSLi). The fine slope is jo which is in agreement with 
our numerical calculations for monolayer films. This suggests that the contribution to the 
eddy currents from precessing magnetization (i.e. from the particular solution of the inho- 
mogeneous Eq.3) is small and the eddy current in the ferromagnetic film is predominantly 
directly induced by the microwave magnetic field of the transducer. The total in the film 
consists then from the transducer Oersted field and the field of the directly excited eddy 
current with negligible contribution from precessing magnetization. 

From the linearized torque equation (8) one finds that the dynamic magnetization is 
driven by the total field hx- If one neglects the small contribution from the precessing 
magnetization to the total field, one can consider (8) as an inhomogeneous equation with 
the right-hand (driving) term in the form of the linear hx{y) function. This inhomogeneous 
system of differential equations can be easily solved analytically, but performing this is out of 
scope of this paper. Here we just mention that if the resonance modes are well-resolved, the 
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resonance amplitude of the i-th resonance should be proportional to the overlap integral 
li = /(f^"*"^^ mW(y) ■ hx{y)dy, where m\^^{y) is the magnetization profile of the i-th eigenmode 
of the film. The latter is calculated as an eigenfunction of the operator which is obtained 
from Eq.(8) by setting a=0, hj.=0. 

Obviously, if the driven resonances are well-resolved (i.e. the resonance linewidths AH ^ 
au/'-f are smaller then distances between the neighboring resonances), the rrix distributions 
in maxima of resonances (Fig 2(a)) are very close to the respective ^«|^(y) and can be used to 
discuss Jj. As one sees from this figure, when the surface current is applied on the Permalloy 
side of the bi-layer film the maximum of the proper distribution of dynamic magnetization 
for the fundamental mode mW(y) coincides with the maximum of the total driving field 
hx{y) and the value of the overlap integral Iq is maximized. On the contrary, if the current 
is applied at the Co-surface of the bi-layer the maximum of the driving field coincides with 
the minimum of the dynamic magnetization. Iq is noticeably smaller, resulting in a much 
smaller amplitude of the fundamental mode. 

This consideration does not apply to the first standing spin wave of the bi-layer. The 
mx{y) distribution for this mode is a quasi-antisymmetric function, but the profile hx{y) 
remains the same as for the fundamental mode (Fig. 2(b)) and is characterized by a large 
anti-symmetric component. As a result the overlap integral Ji does not considerably depend 
on the side at which the microwave current is applied, and its value is large. Thus, the 
behavior of resonance amplitudes in Fig.l is explained not as increase in excitation efficiency 
of the 1st SSW, but as decrease in efficiency of excitation of the fundamental mode for the 
specific bi-layer film orientation with respect to the microwave transducer. 

Since eddy currents are involved in this effect, frequency dependence of the amplitudes 
takes place. Figure 3 shows 5'21/S'21o in the maximum of resonance for different modes as 
a function of frequency. Figure 3(a)) shows the absolute values, and Fig. 3(b) shows the 
relative amplitudes of the standing-wave modes with respect to the fundamental mode. One 
sees that the first exchange mode for Co facing the microstrip becomes dominant at 7 GHz 
or so. 

The model also demonstrates more efficient excitation by the microstrip transducers of 
higher-order modes in conducting monolayer samples than in insulating films (Fig. 4). In 
the case of insulating films, if the surface spins are unpinned d = (Eq.(lO)) the only mode 



which couples to the transducer field is the fundamental mode 



38l] . From Fig. 4 one sees that 
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in the case of conducting films the first SSW which is characterized by an anti-symmetric 
profile m^l may provide considerable response. As seen in our simulation, this response is 
practically frequency independent which is in contrast to the bi-layer films. 

V. CONCLUSION 

In this work we theoretically studied peculiarities of ferromagnetic resonance response 
of conducting magnetic bi-layer films of nanometric thicknesses excited by wide microstrip 
lines. We found strong asymmetry of the response. Depending on ordering of layers with 
respect to the transducer either the first higher-order standing spin wave mode, or the 
fundamental mode showed the largest response. Amplitude of which mode is larger also 
depends on the driving frequency. This theory is in a good agreement with an experiment 
published elsewhere. 

This effect is explained as shielding by eddy currents induced in the film. Our results 
show that for films with thicknesses well below the microwave skin depth this shielding 
remains very efficient. This finding may be useful for studying buried magnetic interfaces, 
as it allows more efficient excitation of higher-order standing spin-wave modes carrying 
information about interface spins. 
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VII. APPENDIX: THE DISCRETE MODEL 

Here we show how we construct the discrete model and how we incorporate boundary 
conditions into it. We demonstrate it using the equation (3) for the dynamic magnetic field 
hx- The discrete version of Eq.(8) with boundary conditions (9) and (10) has a similar form. 
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We use a three-point formula for discrete differentiation to obtain the equation as follows: 

(4^+^) + h^J-'^ - 2h^J^)/A^ + iS-^h^J^ + lATTd-^m^J^ = 0, (22) 

where A is the mesh step along y and j is the number of a point on the mesh. We locate 
the points on the mesh such as no point is at the boundary. In particular the first point on 
the mesh j = 1 is at |/ = A/2, the last point j = n is at y = Li + L2 — A/2. The points at 
the interface of two layers have numbers uq and uq + 1 and are situated at y = Li — A/2 
and y — Li + A/2 respectively. Eq.(23) is valid for any values of j except for 1, n, no, and 
no + 1. For these boundary points boundary conditions should be included into the discrete 
second derivative in Eq.(23). 

We assume that the axis y goes along a horizontal line from the left to the right. Let 
us first consider the point j = Uq to the left of the interface y = Li which is located 
half-way between uq and no -l- 1. We denote an auxiliary point which is located on the 
interface as j = no -|- 0.5. The value of magnetic field at this point /^(."o+o.s) -^^ obtained 
by extrapolating hx{y) dependence beyond the point no from the left using the Taylor 
expansion. First one calculates the field at the point h^°~^'^^ which is half way between no 
and no — 1. Using the Taylor series one obtains: 

h(no-o.5) (23) 
= 4"°~'^ + ('9/ii"°"'V'9l/)(A/2) 
+ {d^h^^°-'ydy^){A/2f 
_ ^(no-i) + (/i(no) _ /i(;^o-2))/4 

+ {h^x°~^^ + 4"°^ - 24"°-i))/8- 
This formula is easily simplified to read 

^(no-0.5) ^ 3^(no) /g + S^"""^) /4 - h^^°-^^ /8. (24) 

Similarly at the point no + 1.5 which is half way between the two first points no -l- 1 no -l- 2 
to the right from the boundary we obtain 

l^{no+1.5) ^ 3/^(no+l) /g + ^4 _ ^(no+3) (25) 

The value of magnetic field at the point at the boundary j = no -l- 0.5 is h^^\ The first 
derivative dh^/dy enters the upper of boundary conditions (4). To evaluate it at the left 
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from the boundary we use the Taylor expansion again 

a/if /dy = dhf^^/dy + {d^h^^'>ydy^)A/2 (26) 

Similarly, to the right from the boundary one has 

dh^^^/dy = (-3/if + 44"°+i) - 4"°+^-5))/A. (27) 

Then substituting Eqs.(27) and (28) into the first of the boundary conditions (4) with Eqs. 
(25) and (26) we obtain 

^(6) ^ (29/i^+^) - 6/i("o+^) + /i^+3))ai ^28) 
^ 24((Ti + (T2) 

(29/i(,"») - 6/i(,"o-^) + /ii"o-2))a2 
24(cTi + aa) ' 

This allows one to evaluate the second derivative at the points no and no + 1. Thus, instead 
of (23) for the point no we have 

A{hf + 4"°-°-^) - 2/ii"°))/^' + iS'^hf'^^ + i47r5-'mi"°^ = 0, (29) 

with /if and /i^'*°~°'^^ defined by Eqs. (29) and (25) respectively. A similar expression is 
easily obtained for no + 1. These finite difference equations now incorporate the inter-layer 
electro-dynamic boundary conditions (4). 

Similarly, at the outer boundary y — Li -\- L2 from (9) we have /if = 0. Then from 
Eq.(30) with uq = n one derives the finite difi^erence equation for j = n 

(13/ii") - 6/if + /ii"-'))/(2A2) + i(^-2/iH + i47r(5-2mf ) = 0, (30) 

where /i^"~°'^) is obtained from Eq.(25) by setting no = n. In the same way for the point 
J = 1 we have 

(-13/i«/2 + 3/ii^) - /if /2)/(2A^) + ir^/if + i47rr^mf = 4r^io, (31) 

where r/ = 1/80 is the factor which relates jo measured in Ampere per meter to measured 
in Oersteds. 

Eq.(8) can be discretized in a similar way incorporating exchange boundary conditions 
(9) and (10) into the exchange operator 2A/M'^^^ at the points j — l,no,no + l,n. In 
particular, for the point n one has 

a^mf ) {dA - 2)(6mf-i) - mf-^)) - {13dA - 10)mf ) 



dy^ 2A2(3 - dA) 
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(32) 



where d is the surface spin pinning parameter from Eq.(lO). 

A — * 

The derived finite-difference equations form a system of hnear algebraic equations Cu — f, 
where u consists of amphtudes of rrix, my, and at the points j = 1, 2...n. The matrix of 
coefficients of this system C has a band form. The system is inhomogeneous, with just one 
element of the vector / being non-zero. This element is the right-hand side of Eq.(32). It 
plays the role of the "excitation term" for the motion of magnetization. The system Cu — f 
can be easily solved using numerical methods of linear algebra. 
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VIII. FIGURE CAPTIONS 



Fig. 1. (Color online) Transmission coefficient 5'21/5'21o for a bi-layer ferromagnetic 
metallic ffim. Layer 1: thickness: 10 nm, saturation magnetization 47rM: 15080 Oe, ex- 
change constant: 1.0 • 10^ erg/cm, Gilbert damping constant: 0.016, conductivity: 1.8 • 10^ 
Sm/m. Layer 2: thickness: 87 nm, saturation magnetization AttM: 8042 Oe, exchange con- 
stant: 0.55 ■ 10^ erg/cm, Gilbert damping constant: 0.008, conductivity: 4.5 ■ 10^ Sm/m. 
Interlayer exchange constant ^412 — 2 • 10^ erg/cm, gyromagnetic constant: 7 = 27r • 2.92 
rad-MHz/Oe, spins at the outer surfaces of the film are unpinned d—0. Microwave trans- 
ducer is a microstrip line 1.5 mm in width, (a) microwave frequency is 4 GHz, (b) 7.5 GHz, 
(c) 18 GHz. Thick lines: amplitudes (left axes); thin hues: phase (right axes). Solid hues: 
Co layer facing the transducer, dashed hues: Py layer facing the transducer. 

Fig. 2 (Color online) Distributions across the film thickness, (a) in-plane component of 
dynamic magnetization rrix- (b) in-plane dynamic magnetic field for the fundamental 
mode, (c) microwave electric field Cz for the fundamental mode. In all panels thick lines are 
for the fundamental mode, thin lines are for the first standing spin wave mode. Black solid 
lines are for the Co layer facing the transducer. Dashed red fine is for the Py layer facing 
the transducer. Blue lines: phase (right axes). All the other parameters are as for Fig. 1. 

Fig. 3. (Color online) (a) absolute values of absorption amplitudes S'21/S'21o as a 
function of frequency, (b) relative amplitudes of standing spin-wave modes with respect 
to the amplitude of the fundamental mode. Bold solid lines: fundamental mode, thin blue 
lines: first standing spin-wave mode. Other thin lines: higher-order standing spin waves 
seen in Fig. 1(c). Solid black and blue lines: Co layer facing the transducer. Dashed red 
lines: Py layer facing the transducer. Parameters of calculation are the same as for Fig. 1. 

Fig. 4. (Color online) Absolute values (right axis, black lines) of absorption amplitudes 
5*21 / 5*210 and their ratio (left axis, dark cyan line) as a function of frequency for a monolayer 
Permalloy film 97 nm in thickness. Thin line: 1st standing spin-wave mode; thick line: 
fundamental mode. All the other parameters are as for the Permalloy layer in Fig. 1. 
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